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Numerical calculation of magnetic form factors of complex shape nano-particles 

coupled with micromagnetic simulations. 

Fatih Zighem 1 ' 2 , Frederic 0tt 2,a [j and Thomas Maurer 3 ' 4 , Gregory 

Chaboussant 2 ' 3 , Jean- Yves Piquemal 5 and Guillaume Viau 6 

1 LSPM, CNRS-Universite Paris XIII, 93430 Villetaneuse, France 

2 CEA, IRAMIS, Laboratoire Leon Brillouin, 91191 Gif sur Yvette, France 

3 CNRS, IRAMIS, Laboratoire Leon Brillouin, 91191 Gif sur Yvette, France 

A LNI0, CNRS-Universite de technologie de Troyes, 10000 Troyes, France 

5 IT0DYS, CNRS-Universite Paris VII, 75 20 5 Paris, France and 
6 LPCN0, CNRS-INSA-Universite de Toulouse, 31077 Toulouse, France 

We investigate the calculation of the magnetic form factors of nano-objects with complex geomet- 
rical shapes and non homogeneous magnetization distributions. We describe a numerical procedure 
which allows to calculate the 3D magnetic form factor of nano-objects from realistic magnetization 
distributions obtained by micromagnetic calculations. This is illustrated in the canonical cases of 
spheres, rods and platelets. This work is a first step towards a 3D vectorial reconstruction of the 
magnetization at the nanometric scale using neutron scattering techniques. 



I. INTRODUCTION 

The recent progress in solid state chemistry has led 
to the possibility of synthesizing nano-objects of non 
spherical shapes: wires by organometallic chemistry [2 J 
or electrochemistry [lj. In particular, the use of the 
polyol process has made it possible to produce well de- 
fined monodisperse magnetic nano-objects. Depending 
on the synthesis conditions, various shapes of particles 
can be obtained such as rods, wires, dumbbells, diabo- 
los, platelets... (see Figure [l]) [3-5J. These are the typical 
forms of magnetic nano-objects which will be the focus 
of this communication. While magnetic nanospheres (as 
found for example in ferrofluids) have been extensively 
studied by Small Angle Neutron Scattering [6HTT]. more 
complex shaped magnetic nano-objects have rarely been 
studied. 

In this communication we focus on the detailed de- 
scription and calculation of neutron scattering on mag- 
netic nano-objects in SANS experiments. We describe 
numerical tools to calculate the magnetic form factors 
of arbitrary shape nano-objects. We present a practi- 
cal procedure which allows to calculate magnetic form 
factors either from an a priori knowledge of the magneti- 
zation distribution or from a minimization of the micro- 
magnetic configuration. 



II. NUCLEAR AND MAGNETIC FORM 
FACTORS CALCULATIONS FOR NEUTRONS 

In this section we first recall the interaction of neutrons 
with small nano-objects. We introduce the nuclear and 
magnetic form factors of nano-objects and describe the 
quantities which can be measured in a Polarized Small 



Angle Neutron Scattering (PSANS) experiment. A prac- 
tical way of calculating the magnetic form factor of com- 
plex nano-objects is then presented. 



A. Nuclear and magnetic form factors 

Let us consider a magnetic nano-object of nuclear Scat- 
tering Length Density (SLD) p (r) which creates an in- 
duction distribution B (r). We define the nuclear form 
factor /n(Q) as the Fourier transform of the nuclear SLD 
distribution, Q = kf — k{ being the scattering wave- 
vector: 



/jv(Q)oc fffp(f)e-^- f df (1) 



Similarly we define the magnetic form factor fs [Q 
as the Fourier transform of the induction field B (r): 



/b(q)« IIJ ' B{r)e- i( l ? dr (2) 

It can be shown that only the component of fs ( Qj 
perpendicular to Q contributes to the scattering |12| : 



Ib± = qx Ib x g = f B - \Jb ■ qj q 



(3) 



with q — jgr being the unit vector along the Q direc- 



tion. 
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The total scattering cross section writes: 
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Figure 1: TEM images of magnetic nano-objects synthesized in various conditions: (a) C080N120 nanowires with rounded tips 
of length 200 nm and diameter 7 nm ; (d) Co^oNi^o dumbbells ; (e) Co diabolos and (f) Co platelets. 



where we recall that a is the neutron spin operator. 

In the case of polarized neutrons with polarization 
analysis, it is possible to measure 4 quantities correspond- 
ing to "Non Spin Flip" (NSF, (++) or (- -)) and "Spin 
Flip" (SF, (+-) or (-+)) scattering as recently experimen- 
tally put into evidence [131 [14] : 



da++ 
dQ 

da~~ 

dn 

da + ~ _ 
dQ 

da~ + 
dQ 



f N (Q) + f B±z ,(Q) 

f N (Q)-fB ± AQ) 

f B±x ,(Q)-ifB ±yf (Q) 

fB ± AQ) + tfB ± AQ) 



(5) 



where f B±a: ,(Q), fB ±y ,(Q) and /b ±z ,{Q) refer to the 
component of fs ± [Qj along the (Ox'), {Oy') and {Oz') 

axis respectively, where {Oz') is the quantification axis of 
the neutron spin defined by the applied magnetic field. 
If the induction distribution is even around the (Ox') 
and {Oy') axes, the Fourier transforms /b ±x ,{Q) and 
Ib ± , (Q) are real so that the two spin-flip cross sections 
are equal. These formulae are derived from Eq. fusing 
the Pauli matrices [15J. 

In the case of polarized neutrons without polarization 
analysis which in practice is the situation encountered in 
most PSANS experiments, the two scattering intensities 
which can be measured are: until now, 2 quantities can 
be measured: 
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(6) 
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In(Q) - f B±z , (Q) + f B±x , (Q) - i f B±y , (Q) 

For non polarized neutron scattering: 



da da~i~ , da 

dQ ~ dQ "•" dQ 



f N (Q)\ 2 + \f B ±(Q)\ 2 (7) 



B. Practical calculation of the nuclear and 
magnetic form factors 

We shall now discuss how to calculate in practice 
the quantity /b±{Q)- A first approach consists in 
assuming an a priori knowledge of the magnetization 
distribution in the particle (for example, M can be 
assumed to be homogeneous). Under this assumption, 
the induction field distribution B (r) can be calculated 
using standard electromagnetic softwares. In the case of 
objects with an axis of revolution, we suggest the use of 
the Femm package [16J. The package provides various 
tools (scripting langage LUA, Octavefemm interface to 
Octave or Matlab, Mathfemm interface to Mathematica) 
which allows extracting the induction field distribution 
~B (r) from the electromagnetic calculation. 

A more general approach, with no assumptions on 
the magnetization distribution in the nanoparticle cal- 
culates the magnetization distribution in the particle un- 
der a given applied magnetic field. Several micromag- 
netic packages have become available for non expert users 
during the last few years {OOMMF [17J, MagPar [18], 
Nmag [19J). We have used the Nmag package since it is 
based on finite elements and is thus especially well suited 
to the types of particles encountered in SANS scatter- 
ing (spheres, cylinders, disks). The procedure consists 
in defining an object geometry and discretization using 
the Netgen mesher [20J (Figure [3]). The 3D magnetic 
moments distribution M (r) in the object can then be 
calculated using the Nmag package (Figure^). 

In order to perform a numerical calculation of the 
magnetic form factors of such magnetic objects, a 
volume of space V containing the objects is defined. Its 
dimensions L x x L y x L z are such that this volume is 
significantly large than the nanoobjects (3-10 times). 
From the magnetization distribution M{r) in the 
nano-object, the induction field B (r) can be calculated 
in the whole volume V. In practice, this induction 
field B {f) is calculated on a regular grid and mapped 
on a 3D matrix of size {n xi n yi n z ). Python scripts 
performing this task are available on the LLB Website 
[21| . This induction distribution is then exported 
into octave or matlab in order to perform numerical 



calculations of the magnetic form factor fs(Q) using 
equations ([2| and |3|. The corresponding scripts are 
freely available at [21J. One obtains a set of three 3D 

matrices {/^(Q), fBy(Q), Ibz(Q)\ describing the 3 
components of the magnetic form factor in the reciprocal 
space. The reciprocal space is mapped with a sampling 
given by AQ = j^ and the useful accessible Q-range 
goes from Qmin = j 1 : to about Q™ ax ^ 0.5 x j^m (e.g. 
for U = lOOnra and ra» = 100, 0.06 < Q < Znm 11 ) 

The Fourier transforms of the induction field are per- 
formed using a FFT algorithm with the condition that 
the dimensions (n x ,n y: n z ) of the matrices are powers 
of 2. Since we are dealing with 3D matrices, the mem- 
ory requirements diverge rather quickly as the matrices 
sizes increase. A matrix size of 128 x 128 x 128 provides 
both a reasonable memory footprint (32 MBytes / ma- 
trix) and fast calculations times (507718 per FFT on a 
standard desktop computer). Such a sampling rate pro- 
vides accurate Fourier transforms in the reciprocal space 
over 2 orders of magnitude (see Figure 5 for example). 
We found that a volume of space V about 6 times larger 
(i.e. L x = 6D for a sphere of diameter D) than the 
studied object provided good results in the Q range of 
interest. 

We now define the scattering geometries that will be 
considered in the actual numerical calculations presented 
in the next section. We first underline that in the general 
case of polarized neutron scattering on anisotropic parti- 
cles, it is necessary to define 3 axis systems. The first one 
(x, y, z) is attached to the studied nano-objects. The sec- 
ond one (X, Y, Z) describes the spectrometer geometry. 
The third one (x',y',z f ) describes the orientation of the 
quantification axis (defined by the applied magnetic field 
along z'). In order to describe a scattering experiment, 
it is necessary to define the relative orientations of these 
3 sets of axis. 

In order to simplify the presentation and to com- 
ply with experimental conditions, we consider a scat- 
tering plane (Qx = Qcosa, Qy = Qsina) correspond- 
ing to the 2D detector plane of a SANS spectrometer 
(XOY). For each scattering direction a in the scatter- 
ing plane, the perpendicular components of the magnetic 
form factor f± (0, r) can be calculated using equation |3|: 
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for a; = 0, a cylinder is aligned along the (OZ) axis. 
For uj = 90°, a cylinder is aligned along the (OX) axis. 
The distribution B (f) XYZ in the spectrometer axis is 
obtained by a rotation uj : 



(9) 



For the sake of illustration, we shall not include both 
the nuclear scattering and the magnetic scattering in the 
calculations. Even though this is not realistic in a gen- 
eral scattering experiment where the nuclear contribu- 
tion is usually large, such situations can be encountered 
in SANS experiments: the nuclear contrast of cobalt 
particles can easily be matched by using an appropri- 
ate deuterated solvent [TJ |H] ; ferromagnetic particles in a 
paramagnetic matrix also allow to extinguish the nuclear 
contrast [TT]. Ignoring the nuclear scattering effects will 
allow us to focus on the magnetic effects which are the 
scope of this communication. 

In the following, we shall different relative orientations 
of the applied magnetic field direction with respect to the 
scattering plane: 

1. Longitudinal scattering (H || OZ) . Equation ([5| 
becomes: 



da ++ _ da_ 

dn 



dQ 



da^ 



da 



dn 



dn 



fB ±z (Q) 



fB ±x (Q)TifB ± AQ) 



(10) 



2. Transverse scattering (H \\ OX) . Equation ^ 
becomes: 



da++ 
dQ 

da + ~ 
dQ 



da 
dQ 

da~ + 
dQ 



f Bxx (Q) 



fB ±Y (Q)TifB ±z (Q) 



(11) 



Note that even though the above formulae are very simi- 
lar, the key parameter is the distribution fs ± (Q) with re- 
spect to the scattering plane. The scattering geometries 
(l)-(4) presented on Figure [2] will lead to qualitatively 
different results which are discussed in the next section. 




fs x ~ (fs x cos a + Ib y sm a ) cos a 

Ib y ~ Ub x cos a + fs Y sin a) sin a 

Ib z 

(8) 

Note that the above formula is only valid in the 
{Qx, Qy) scattering plane. In the practical procedure, 
the magnetic induction distribution is calculated in the 
(xyz) reference frame where (Oz) is the revolution axis 
of the object. The examples proposed in [21 J follow this 
convention. We define uj as the rotation angle of the 
object around the (OY) axis. With these conventions, 



III. MODEL CALCULATIONS 

In this section we will present the calculation of the 
magnetic form factors of different types of nano-objects: 
(a) a sphere of diameter D = 10 nm, (b) a flat disk of 
diameter D = 10 nm and height h = 2.5 nm and (c) a 
cylinder of diameter D = 10 nm and length L = 50 nm. 
These objects are presented in Figure [3] The magnetic 
parameters used are the following: magnetization satura- 
tion of cobalt M = 1.4 x 10 6 A.m _1 (=1.76 T); exchange 
stiffness A = 1.2 x 10 11 J.m -1 . Moreover, a uniaxial 



Longitudinal scattering 



Transverse scattering 




(2) 



ESfe 



XT' 




Figure 2: The different scattering that will be considered. The neutron beam is incident along the (OZ) direction. (XY) defines 
the scattering plane. The longitudinal (1-2) and transverse (3-4) scattering geometries correspond to a magnetic field applied 
along (OZ) and (OX) respectively. The scattering object can be oriented with its easy axis parallel (1-4) or perpendicular 
(2-3) to the applied field. 




Figure 3: Typical meshes obtained from the Netgen mesher for the different nano-objects: (a) a sphere of diameter D — 10 
nm, (b) a flat disk of diameter D = 20 nm and height h = 2.5 nm, (c) a cylinder of diameter D — 10 nm and length L — 50 
nm, (d) a disk of diameter D — 50 nm and height h — 5 nm. Note that the distance between two nodes is around 1 nm. (e) 
Maps of the magnetic moments distribution at remanence for the disk defined in (d). The colors encode the ^/-component of 
the magnetization (m y — M y /M) in the profile view (left) and z-component of the magnetization (m z — M z /M) viewed from 
the top (right). 



magnetocrystalline anisotropy K = 7 x 10 5 J.m -3 along 
z was considered in the case of the flat disk (Figure [3b ) 
in order to force the magnetic moments to be normal 
to the disk surface. In a first step, the remanent con- 
figuration was determined for all the objects. It reveals 
a rather homogeneous distribution of the magnetic mo- 
ments along the Oz direction for objects (a), (b) and 
(c). This is due to the fact that these objects are to the 
first order approximations of elliptical objects in which 
the demagnetizing field in homogeneous. For simplicity, 
only nano-objects with a quasi uniform magnetic mo- 
ments distribution have been considered. However, more 
complex distributions can be considered as it is highlight 
for a disk of diameter D = 50 nm and height h = 5 nm 
within which a vortex state appears (see Figure^). 

We perform numerical calculation of the 3D nuclear 
form factors so as to be able to compare the results 
with analytical formulae and validate the numerical pro- 
cedure. As a next step, we perform calculations of the 



magnetic form factors where we illustrate the contribu- 
tions of the demagnetizing and the dipolar fields. 



A. Validation of the numerical procedure: 
application to the calculation of nuclear form factors 

Although this paper is devoted to the magnetic form 
factors of nano-objects, we have also numerically calcu- 
lated the 3D nuclear form factor of these objects in or- 
der to check the applicability of our numerical algorithm. 
When calculating the 3D nuclear form factor, the volume 
of space V containing the object is mapped onto a 3D 
matrix G x , y , z describing the object geometry in such a 
way that G x ,y,z — 1 fc> r (%, V, z) inside the object and 
Gx,y,z = for (x, y, z) outside the object. 

Figures [4^i-b-c present the 3D nuclear form factors cal- 
culated for the sphere (D = 10 nm), the flat disk (D = 10 
nm, h = 2.5nm) and the cylinder (D = 10 nm, L = 50 
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Figure 4: Numerically calculated form factors (circles) and comparison with analytical formulas (continuous lines), a), b) and 
c) correspond to 3D form factors calculated for the sphere D — 10, the cylinder (D — 10 nm, L — 50 nm), the flat disk (D — 10 
nm, h — 2.5 nm) and ID cuts of the form factor in the reciprocal space in specific directions. 



nm) defined in Figure [3] The scattering is isotropic in the 
sphere case and anisotropic in the other cases. The nu- 
clear form factors of these simple objects are well known 
and can be analytically calculated [22J. We can thus 
compare our numerical results with the analytical formu- 
lae. Figures [4]i-e-f compare the numerical and analytical 
calculations along the main directions in the reciprocal 
space (in nm -1 ) for the different objects. The agreement 
is satisfactory for all the objects which proves the appli- 
cability of our numerical calculations. 

These results show that it is possible to calculate rather 
accurately, in a single shot, the form factor in a Q-range 
extending over 2 decades with scattering intensities ex- 
tending over 3 decades. If one would be interested in 
a different Q-range, it is necessary to resize the V box 
around the object. As Qmin = j^, m order to probe 
smaller Q values, the size Li of the box should be in- 
creased As Q™ ax no 0.5 x j^rii, in order to probe larger 
Q values, either the size Li of the box should be decreased 
of the number of mapping points rii should be increased. 
This is similar to a real SANS experiment where the mea- 
suring conditions are changed in order to cover a wider 
Q-range. 

In the following we shall however restrict ourselves to 
a V box 6x6x6 times the size of the object, which 
provides an appropriate sampling in the Q space for our 
demonstrations. 



B. Case of magnetic objects 

The same calculation procedure was applied to the case 
of magnetic objects. Figure 5 presents the calculation of 
the magnetic form factors in the case of the 3 canonical 



geometries (sphere, disk and cylinder). A first reference 
calculation was performed by calculating the form factor 
of an homogeneously magnetized sample and for which 
only the magnetization inside the sample was considered 
(blue lines on Figure 5). A second calculation was per- 
formed where the form factor of the induction map was 
also included. The induction map is plotted on Figure 
5 a-c. We considered the induction in the sample given 

by B int .(r) = Mo \M(r) + H d (r)j where H d (r) is the de- 
magnetizing fields and the stray fields outside the sample 
(B ext \r)). The contribution of these inductions distribu- 
tions was calculated separately (black and green curves). 
The sum of these contributions B tot \r) overlaps with the 
form factor simply obtained from the magnetization dis- 
tribution (see figure 5d-e-f). The deviation observed at 
low Q are numerical artifacts due to the fact that the 
stray fields expand beyond the calculation box. This 
numerically demonstrates that the magnetic form factor 
can be obtained by considering the magnetization distri- 
bution without taking into account the demagnetization 
fields and the stray fields outside the sample. 



IV. CONCLUSION 

We have developed numerical tools [21 J which per- 
mits the calculations of the structural and magnetic form 
factor of nano-objects with complex geometrical shapes. 
Such calculations were until recently limited due to the 
extensive memory usage of 3 dimensional FFT calcula- 
tions. As Polarized SANS spectrometers are becoming 
widely available the problem of quantitatively process- 
ing the data obtained on complex magnetic nano-systems 
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Figure 5: Longitudinal scattering on 3 types of objects (sphere, cylinder and disk) magnetized along the applied field, a), b), 
and c) Induction distribution in (Y Z) plane for a sphere, a cylinder and a disk magnetized along the Z direction, d), e) and f) 
Scattering cross sections of different magnetic contributions: the magnetization only (M(r)) } the internal magnetic induction 
(Bi n t.(r))i the dipolar field created by the magnetic moments outside the object (B ex t.( r )) an d the whole magnetic induction 
(B tot .(f) = Bint.fi) + Bext.i)). 



will arise. This first step will allow experimentalists to 
compare their measurements with realistic models of the 
magnetization in their nano-systems. In the future, with 
the widespread availability of systems with tens of gi- 
gabytes of memory, one may consider applying similar 
methods for the numerical calculation of complex struc- 
ture factors. 
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